What are universal features of gravitating Q-balls? 
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We investigate how gravity affects Q-balls by exemplifying the case of the Affleck-Dine potential 
V{4>) :— m ln(l + ^j). Surprisingly, stable Q-balls with arbitrarily small charge exist, no matter 
how weak gravity is, contrary to the case of flat spacetime. We also show analytically that this 
feature holds true for general models as long as the leading order term of the potential is a positive 
mass term in its Maclaurin series. 
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I. INTRODUCTION 

Q-balls a kind of nontopological solitons [2J, ap- 
pear in a large family of field theories with global U(l) 
(or more) symmetry. It has been argued that Q-balls 
with the Afflcck-Dinc (AD) potential could play impor- 
tant roles in cosmology [3j. For example, Q-balls can be 
produced efficiently and could be responsible for baryon 
asymmetry Q and dark matter Therefore, their sta- 
bility is an important subject to be studied. Because 
Q-balls are typically supposed to be microscopic objects, 
their self-gravity is usually ignored; and accordingly, sta- 
bility of Q-balls with various potentials has been inten- 
sively studied in flat spacetime [6|-|9|. 

However, if we contemplate the results on boson 
stars [To|, we notice the possibility that self-gravity can 
be important even if it is very weak. For example, in the 
case of the potential Vmini = m 2 (/> 2 /2, equilibrium solu- 
tions, called (mini-)boson stars, exist due to self-gravity, 
though no equilibrium solution exists without gravity. 
An important point is that there is no minimum charge 
for (mini-)boson stars, i.e., they can exist even if their 
self-gravity is very weak. This example suggests the im- 
portance of the unified picture of Q-balls and boson stars, 
which are different from each other solely in model pa- 
rameters. 

From this motivation, we have investigated the follow- 
ing three models. The first one is [ll[ 



W) := 



H(f> 3 + \(f) 4 with m 2 , A>0, (1.1) 



which describe Q-balls (// > 0) and boson stars (/i = 0) 
comprehensively. The second one is [l2j 



Vi(4>) := ^-0 2 -A0 4 + vW with m 2 , A, M > 0. (1.2) 
2 M z 
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Interestingly, gravitating Q-balls (boson stars) in these 
two models have a common property that stable Q-balls 
with arbitrarily small charge exist no matter how weak 
self-gravity is, while Q-ball properties in flat spacetime 
fairly depend on the potentials. 

Although these examples reveal gravitational effects 
on Q-balls, both models are described by polynomial in 
<f>. From the theoretical point of view, however, models 
with logarithmic terms are more natural. Specifically, in 
the AD mechanism noted above, there are two types of 
potentials: gravity-mediation type and gauge-mediation 
type. The former type is described by the potential, 



^grav.(0) 



1 + Jfln 



M 



with 



M > . 



(1.3) 



In flat spacetime equilibrium solutions for K > are 
nonexistent while those for K < are existent. If we 
take self-gravity into account, stable Q-balls exist even 
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our previous paper [l<jj |. we have shown that gravitating 
"Q-balls" exist, which are surrounded by Q-matter, even 
for K > 0. Unfortunately, since we cannot define the 
Q-ball charge in this case, it is difficult to say about the 
common property noted above and their stability has not 
been explored yet. 

In the present paper we extend our analysis to the 
gauge-mediation type, 



V'gaugc^) := m 4 In ( 1 + ) with m > . (1.4) 
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We shall show that gravitating Q-balls with (|1.4j) have a 
common property with U m ; n i, V3 and V4, and discuss the 
reason. 

This paper is organized as follows. In Sec. II, we de- 
rive equilibrium field equations. In Sec. Ill, we show 
numerical results of equilibrium Q-balls. In Sec. IV, we 
explain why self-gravity gives a rather model indepen- 
dent feature even if it is very weak. In Sec. V, we devote 
to concluding remarks. 
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II. ANALYSIS METHOD OF EQUILIBRIUM 
Q-BALLS 



A. Equilibrium field equations 

We begin with the action 
S = / d A x^J —g 



16ttG 



where <p 



(2.1) 

62) is an S0(2)-symmetric scalar held 



and (f> := \J 4> ■ d> = We assume a spherically 

symmetric and static spacetimc, 



ds 2 



-a 2 {r)dt 2 +A 2 (r)dr 2 +r 2 (d9 2 +sm 2 ddip 2 ). (2.2) 



For the scalar field, we assume that it has a spherically 
symmetric and stationary form, 



(01) 4>2) = </>(r)(cos wi, sinwi). 
Then the field equations become 
rA 3 



(2.3) 



-G := 



a'+^V-i) 



0^ , w 2 2 

2A 2 



2a 2 



V 



VOL 
2 ^rr 



= 47rGraA 



0' 2 w 2 2 



.4 2 



-□c 



2A 2 2a 2 
2 a' A' 



r 



.4 



ujA 
a 



(2.4) 



(2.5) 



(2.6) 



where ' := d/dr. To obtain Q-ball solutions in curved 
spacetime, we should solve (|2.4j) - (l2.6[) with boundary 
conditions, 

A(0) = A(oo) = a(oo) = 1, 

A'(0) = a'(0) = ^(0) = <?!>(oo) = 0. (2.7) 

We also restrict our solutions to monotonically decreas- 
ing <p(r). Because of the symmetry, there is a conserved 
charge called Q-ball charge, 



d c xi/^gi f v {(j}id v <p2 - <p2d v <pi) = ul, 



, Ar 2 ^ 2 , 
where I := An I dr. 



(2- 



We suppose V'gauge Model \1A\ . Rcscaling the quanti- 
ties as 



= -, v g 



^gaugo 



m 
to 



augc 



^=M1 + ^), 



uj = — . k = Gm , t = mt, f = mr, (2.9) 
m 



the field equations (|2.4I) - (|2.6I) are rewritten as 



A' + -(A 2 1) = AnKrA + + U 



2A 2 • 2a 2 ' ' 
(2.10) 



a' + |(1 - A 2 ) = An.raA 2 ( IL + ^ - W ) , 



(2.11) 



+ 



2 o!__X_ 
f a A 



&A\ ~ _ 2 rfVgaugc 



(2.12) 



B. Equilibrium solutions in flat spacetime 

In preparation for discussing gravitating Q-balls, we 
review their equilibrium solutions in flat spacetime (k = 
0). The scalar field equation (|2.12l) reduces to 



6" = -?0'-c^+^fel 



(2.13) 



This is equivalent to the field equation for a single 
static scalar field with the potential V w := Vg aU g C — 
w 2 </> 2 /2. Equilibrium solutions satisfying boundary con- 
ditions ()2.7[) exist if 

min(K,) < ^gauge(O) and ^(0) > 0. (2.14) 



We obtain 



d 2 V^ 



dV u 20 
di 



<j) i + 4> 2 
2 



= ui 2 c 



(1 + 2 ) 2 1 + . 



(2.15) 
(2.16) 



The first condition in f|2 . 14[) is trivially satisfied since 
is unbounded from below. If we introduce e 2 := 2 — ui 2 , 
the second condition in (12. 14)) leads to 



e 2 > 



(2.17) 



The two limits e 2 — > 2 and e 2 — > correspond to the 
thin- wall limit and the thick- wall limit, respectively. 

If one regards the radius r as "time" and the scalar 
amplitude <p(r) as "the position of a particle" , one can 
understand Q-ball solutions in words of Newtonian me- 
chanics. Equation (|2.13l) describes a one-dimensional 
motion of a particle under the conserved force due to 
the potential — VL(0) and the "time" -dependent friction 
-{2/r)d<j)/dr. 

To discuss gravitational effects later, it is useful to es- 
timate the central value 4>o := 0(0) m fl&t spacetime. 
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Because V u « at r = 0, its order of magnitude is es- 
timated as a solution of V u = (0(0) 7^ 0). For V gauge 
with the thick- wall condition e 2 <C 1, we obtain 



In this case, V u and e 2 should be redefined by 



1 

2 



(2.18) 



where we have used Maclaurin expansion of ln(l+0 2 ) and 
neglected higher order terms O(0q). Then, we obtain 

O ~ e. (2.19) 

which verifies our assumption that higher order of 0o is 
negligible. 



III. GRAVITATING Q-BALLS 
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FIG. 1: For a thick-wall Q-ball for Cj 2 ~ 1.999, we show 
(a) behavior of the metric functions for a gravitating one for 
k = 0.0005 and (b) comparison of 4>(r) with that for the flat 
case K = 0. 

The potential picture described above is effective in 
arguing equilibrium solutions also in curved spacetime. 



V„ := V„ 



augc 



2a 2 



2-^. (3-D 



a* 



r.2 



"The potential of a particle", — K,, is now "time"- 
dependent, which sometimes plays an important role, as 
we see below. 

First, we show the result on a thick-wall Q-ball with 
1.999 in Fig. CD We put k = 0.0005 for a gravi- 
tating Q-ball. The metric functions a(r), A(r) and the 
field amplitude 0(f) are shown in (a) and (b), respec- 
tively. Because a(f) and A(f) are close to one, one may 
think that gravity acts as small perturbations. However, 
looking at 0(f), we find that gravity changes its shape 
drastically. Near the origin, we find that the scalar field 
for the gravitating case takes much smaller value than 
that for the flat case. 




FIG. 2: -Vu, for a thick-wall Q-ball u 2 ~ 1.999. For the fiat 
case k = 0, — V w is fixed as shown by a solid line. On the 
other hand, for the gravitating case k = 0.0005, — Vu changes 
from a dotted line to a solid line as "time" f goes. As a 
result, the scalar field should start from <j> ~ 9.68 x 10 -2 and 
4> ~ 4.96 x 10 -3 to satisfy <f>{po) — for the flat case and the 
gravitating case, respectively as shown in this figure. 



We explain the reason for this result by using the po- 
tential — Vu, as shown in Fig. [3J For the flat case k = 0, 
— Vu, is fixed as shown by a solid line. As a result, the 
scalar field with a relatively large value (cf> ~ 9.68 x 10~ 2 ) 
at the initial time r = rolls down and climb up the po- 
tential and finally reaches = at the time f — > 00. On 
the other hand, for the gravitating case n = 0.0005, —V u 
changes from a dotted line to a solid line as "time" f goes. 
An important point is that the sign of — changes near 

the origin. Therefore, the scalar field should take a suffi- 
cient small value (0 ~ 4.96 x 10~ 3 ) to satisfy 0(oo) = 0, 
as shown in Fig. [5] 

Second, we show the result on a relatively thin Q-ball 
with uj 2 ~ 0.35 in Fig.[3j The metric functions a(r), A(f) 
and the field amplitude 0(f) are shown in (a) and (b), 
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FIG. 3: For a Q-ball for Cj 2 ~ 0.35, we show (a) behavior of 
the metric functions for a gravitating one for k = 0.0005 and 
(b) comparison of <f>(r) with that for the fiat case k = 0. 



respectively, (a) tells us that, compared with the case 
Cj 2 ~ 1.999, the gravitational field is fairly stronger near 
the origin but approaches the flat spacetime faster as 
f — » oo. (b) indicates that the size of the Q-ball becomes 
smaller due to self-gravity. 

We can also explain the reason for this result by using 
the potential — V Ul as shown in Fig. |4j For the flat case 
K = 0, — K; is fixed as shown by a solid line while — for 
a gravitating case k = 0.0005 changes from a dotted line 
to a solid line as "time" f goes. An essential difference 
from the former thick- wall case is that this change occurs 
"quickly" . As a result, the scalar field dp ~ 10.68 (~ 9.85) 
at r = rolls down the potential and finally reaches = 
at f — y oo for the flat (gravitating) case. Thus, (j> for the 
gravitating case rolls down the potential faster than that 
for the flat case. 

Larger influence caused by weak gravity for the thick- 
wall case might seem paradoxical. To understand it, it is 
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FIG. 4: -V u for a Q-ball cj 2 ~ 0.35. For the fiat case k = 0, 
— Vui is fixed as shown by a solid line while — VL for a gravitat- 
ing case k — 0.0005 changes from a dotted line to a solid line 
as "time" r goes. Since this change occurs "quickly", their 
difference is small comparing with that might be expected. 
I.e., the scalar field dp ~ 10.68 (~ 9.85) at f = rolls down 
the potential and finally reaches <j> — at f — > oo for the flat 
(gravitating) case. 



convenient to define the density p and the radial pressure 
p r of the scalar field in the fluid form as 



2A2 + 2a 2 



•"gauge j 



Pr ■ = 



2A 2 2a 2 



- V, 



gauge ■ 



(3.2) 



(3.3) 



Then, let us remember hydrostatic equilibrium equa- 
tions, which are another expression of Einstein equations, 
as in the usual star. 



dm(r) , 
= 4irr 2 p, 



dr 



dp r 
dr 



Gm(r)p 



(3.4) 



(3.5) 



where m(r) is the mass function of the Q-ball. We should 
notice that the pressure gradient must work as a repulsive 
force against the gravity to support the Q-ball while it 
works as an attractive force in the flat case. If we pay 
attention to the values of r and dp in Figs. Q] and [3] (b), 
we find that the absolute value of the pressure gradient 
for the thick-wall case is far smaller than that for the 
thin- wall case. Thus, weaker gravity does not necessarily 
mean a smaller influence to Q-balls and it is interesting 
to investigate its influences for various Cj 2 which will be 
discussed below. 

As we discussed in our previous papers , stabil- 

ity of Q-balls can be easily understood from the relation 
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FIG. 5: (a) Q-E relation for k = 0, 0.0005 and 0.004 and (b) 
its magnification around Q ~ fOO. 
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FIG. 6: Q-(E — |Q) relation corresponding to Fig. 0(b). 



1 


D 


i 


1 

c~ 


fe B 






- = 


: \ 
a\ 






_ 










1 kfO.004 








/k=0.0005 


i i 


.....K=0 

l 





200 400 600 800 1000 1200 

Q 



FIG. 7: Q-<? relation for n = 0, 0.0005 and 0.004. 



between Q and the Hamiltonian energy E, which is de- 
fined by 



r^oo 2GA 



Ms 
2 ' 



(3.6) 



where Ms is the Schwarzschild mass. Here, stability 
means local stability, that is, stability against small per- 
turbations. We also normalize E as 



E := 



E 



(3.7) 



We show Q-E relation for k = 0, 0.0005 and 0.004 in 
Fig. [5] (a) and its magnification around Q ~ 100 in (b). 
Since complicated structures are concentrated near the 
line E = |Q in (b), we also show the corresponding Q- 
(E — jQ) relation in Fig. [SJ We find that this relation 
is similar to that for the V4 model (see, Fig. 3 in 12]). 
Therefore, the stability for the present model can be also 
understood in the same way as that for the V4 model. 

For the fiat case represented by the dotted line, there 
are double values of E for a given value of Q. These two 
branches merge at the point B (Q ~ 77) as shown in 
(b) or Fig. [B] We can understand that the upper(lower) 
branch represents unstable (stable) solutions. 

Next, we discuss the case of k = 0.0005. Figures[5]and 
Fig. [S] tells us that the solution sequence A-B-C is anal- 
ogous to the sequence in the flat case. We can therefore 
interpret the upper (lower) branch from A to B (from 
B to C) represents unstable (stable) solutions. On the 
other hand, we can find the intrinsic differences from the 
flat case from A to the origin and from C to D (and 
sequences of cusp structures). By energetic (or catas- 
trophic) argument described in our previous paper [l2j], 
the branch from A to the origin is regarded as stable. 
The point C corresponds to the maximum of Q (Q max ). 
This suggests that a Q-ball with larger charge than Q ma x 
cannot support itself due to the self-gravity. Both this 
interpretation and energetic (or catastrophic) argument 
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indicates that the solutions C to D (and sequences of 
cusp structures) are unstable. 

As for the case K = 0.004, Q ma x becomes very small 
due to the large self-gravity and there is no fine structure 
like the sequence A-B for k = 0.0005. The smooth curve 
from the origin to the point C represents stable solutions. 

As we discussed in (TT| - [l3| . cl> 2 or e 2 is a state vari- 
able, while Q and k are control parameters, in words of 
catastrophe theory. Therefore it is instructive to depict 
the Q-e 2 , too, in Fig. [7] For the thick- wall solutions 
e 2 <C 1 with gravity, as we explained using Fig. [2j Q-ball 
charge is very small because 0(0) <C 1, and there is no 
lower bound of Q, contrary to the flat case. The solution 
near the point C was also explained using Fig. 0] In this 
case, the Q-ball charge becomes small due to the strong 
gravity. These phenomena occur for \g rr — 1| ~ 1 at its 
peak [ll|, Ojj]. For k = 0.004, the thick- wall regime and 
the strong gravity regime overlap each other. Therefore, 
the solution sequence is quite different from that for the 
flat case. 




FIG. 8: Values of Q for various n corresponding to the points 
A, B, C and D in Fig. [5] 

Because the extremal points A, B, C and D in Fig. [5] 
indicate the main feature of the solutions sequences, we 
also depict how the values of Q at these points vary with 
k in Fig. [SI The values of Q at C, which corresponds 
to Qmax, become smaller as k becomes larger; this is a 
common feature with V3 and V4 models [TTI.Il^. Hij . The 
values of Q at the points A and B merge at k ~ 0.0035 
and the small structure disappears for k > 0.0035, as 
already seen for k — 0.004 in Fig. [5] These continuous 
change as k is easily understood. 

This figure, however, reveals a nontrivial result. The 
local maximum A does not disappear for k — > +0 while 
it is nonexistent for k = 0. This means that there is no 
lower bound of Q for k — > +0 while there is a minimum 
of Q for k = 0. This result is against our naive idea that 
gravitational effects should vanish in the limit of n — > 
+0. Therefore, It should be argued carefully whether or 



not there remain thick-wall solutions with Q —> in the 
weak gravity limit k ^ 0. This is the subject of the next 
section. 



IV. THICK- WALL SOLUTIONS FOR k / 

We consider thick-wall solutions (e 2 <§; 1) with weak 
gravity by expressing the metric functions as 

a 2 = l+h(r), A 2 = l+/(r), (\h\ « 1, l/l « 1), (4.1) 

and we shall take up to first order in h and / hereafter. 
As we discussed for the flat spacetime in Sec. IIB, we 
evaluate 0o := <fi(0) as a solution of V u — with (|4.ip . 



+ £o 2 h{0)-<p 2 ~0, 



(4.2) 



where we have neglected higher order terms 0(4>q). 

Let us consider the limit e 2 — > 0. For the flat case n — 
0, <fio ~ e since h(0) can be taken to be zero identically. 
However, for any small k ^ 0, it is not evident whether 
or not h(0) can be negligible and we should compare the 
order of h(0) with that of e 2 . For this purpose, we should 
also estimate it by using the Einstein equations. To do 
this, we assume a top-hat configuration, 

C 

4>{r) ~ 0o < 1 for r<— , C = const. (4.3) 

e 

For the flat case, 4>q — e, as in Eq. (|2.19p . 
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FIG. 9: e 2 -0o relation for various k. 



From the Einstein equations, we find 



G\ + G\ :-- 



A 



1 . ( w 2 2 \ 
8irKf 2 Aa — - V , 



where i denotes spatial components. If we take the weak 
field approximation (|4.1[) and the thick-wall approxima- 
tion e 2 -C 1, we obtain 



(f 2 h')' ~ Wirnf 2 ^ 2 



(4.5) 
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With the central boundary condition h'(0) = and the 
approximation (|4.3p . we can integrate (|4.5|) as 



ti 



16 

y 



TTK(f) f, for f < 



C 
e 



(4.6) 



With the outer boundary condition h(^f) ~ h(oo) = 0, 
we can integrate (|4.6|) as 



h(0) 



From (14.21) and (14.71). we obtain 



3c 4 



,C 2 



87tkC 2 + 3e 2 ' 



(4.7) 



(4.8) 



This formula clearly shows how 0o converges as e — > 
and k — > 0; it depends on their convergent rates. If 
e 2 ^> kC 2 , we have <fio ~ e, as in the flat case. On the 
other hand, if e 2 -C «C 2 , we have 



2C* V 2ttk 



(4.9) 



In the real situation, k is very small but a nonzero con- 
stant, while e is variable and determined by initial con- 
ditions. Therefore, if we discuss the thick-wall limit in 
weak gravity, the latter result (|4.9|) applies. In this case, 
because 



2^3 



Q oc 4> t 



(4.10) 



there is no lower bound of Q as expected. Moreover, the 
point A in Fig. [7] can be interpreted as the point when 
e 2 ~ kC 2 . To confirm the above argument, we show the 
numerical relation e 2 -(j>o in Fig. |H1 We find that (I4.9[) 
holds true and C ~ O(10). 
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FIG. 10: e 2 -(— h(0)) relation for various k. 

To check the consistency with the assumption of weak 
gravity, we also estimate /i(0) by substituting (14. 9p into 



ft(0) ~ -e 2 . (4.11) 
Because of the thick- wall assumption e <C 1, it is consis- 
tent with the first assumption \h\ <C 1, and confirmed by 
the numerical relation e 2 -(— h(0)) in Fig. [lOj This result 
tells us that h(0) cannot be negligible in Eq. (|4.2p inde- 
pendent of k 7^ 0. Thus, we should distinguish between 
k = (i.e., h(0) = 0) and any small value of k ^ 
clearly. In summary, if e 2 -C kC 2 and e<l arc satisfied, 
we conclude that effects of self-gravity cannot be ignored 
no matter how weak the gravity is. 

We want to examine what types of potentials we can 
apply the above results. To obtain (|4.8[) . we have as- 
sumed that the second leading order is quartic, </> 4 . How- 
ever, to obtain (|4.9[) . which is the result for e 2 <C kC 2 , 
we have only assumed that the leading order in V u (or 
that in the Maclaurin series) is quadratic, 4> 2 . Therefore, 
the above results hold true not only for the model (|1.4j) 
but also general models in which a positive mass term is 
a leading order one. 



V. CONCLUSION AND DISCUSSION 

We have investigated how gravity affects Q-balls by 
exemplifying the case with the AD potential V((f>) := 
m 4 ln(l + ^-f). Surprisingly, stable Q-balls with arbitrar- 
ily small charge exist, no matter how weak gravity is, 
contrary to the case of flat spacetime. The result for the 
gravitating case is a universal property which has been 
known to hold for various potentials |T0l - ll2l |. We have 
also showed that this feature of gravitating Q-balls holds 
true for general models as long as the leading order term 
of the potential (or that in the Maclaurin series) is a pos- 
itive mass term. Therefore, this result does not change if 
we include nonrenormalization terms, which we have ig- 
nored in the gauge-mediated AD potential. Our results 
suggest that gravity may play an important role in the 
Q-ball formation process. 
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